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^ ' Abstract 

X5 : 

QJ We present a sufficient condition for approximate controllability of the bilinear discrete- 

spectrum Schrodinger equation exploiting the use of several controls. The controllability 
OO . result extends to simultaneous controllability, approximate controllability in H'^, and 

tracking in modulus. The result is more general than those present in the literature even 
in the case of one control and permits to treat situations in which the spectrum of the 
U: uncontrolled operator is very degenerate (e.g. it has multiple eigenvalues or equal gaps 

i among different pairs of eigenvalues). We apply the general result to a rotating polar lin- 

r~| ' ear molecule, driven by three orthogonal external fields. A remarkable property of this 

■ model is the presence of infinitely many degeneracies and resonances in the spectrum 
preventing the application of the results in the literature. 
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^ ; 1 Introduction 



(N . 

Q . In this paper we study the controllability and the tracking problem for a multi-input bilinear 
^ I Schrodinger equation 

= {Ho + u,{t)H^ + ... + Up{t)Hp)ij{t) (1) 

[ where Hq, . . . , Hp are self-adjoint operators on a Hilbert space "H and the drift Schrodinger op- 
erator Ho (the internal Hamiltonian) has discrete spectrum. The control functions Mi(-), . . . , Up{ 
representing external fields, are real- valued and ■?/'(■) takes values in the unit sphere of "H. 
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The controllability of system (JT]) is a well-established topic when the state space "H is finite- 
dimensional (see for instance |D'A08) and reference therein), thanks to general controllability 
methods for left-invariant control systems on compact Lie groups ( |Bro72| IJS72| [JK8H IGB82| 
IMGK96]). 

When % is infinite-dimensional, it is known that the bilinear Schrodinger equation is not 
exactly controllable (see |BMS82l [TurOO] ). Hence, one has to look for weaker controllability 
properties as, for instance, approximate controllability or controllability between eigenstates 
of the internal Hamiltonian Hq (which are the most relevant physical states). In certain 
cases, when the space variable is one-dimensional, a description of reachable sets has been 
provided (see |BC06t lBL10| ). In dimension larger than one or for more general situations, the 
exact description of the reachable set appears to be more difficult and at the moment only 
approximate controllability results are available. Most of them are for the single-input case 
(see, in particular, |CMSB09l IMirOQl INerOQl INerlOl [BNTOl IBCCS121 INN12| ). Such results are 
based on sufficient conditions for controllability that are generic |PS10| IMSlOj INerlO] even in 
the case p = 1. Nevertheless, in many examples interesting for applications these conditions 
cannot be directly applied or controllability fails to hold, as a consequence of the symmetries 
of the system. Symmetries can induce degeneracies in the spectrum (e.g. multiple eigenvalues 
or presence of equal spectral gaps) and reduce the coupling of eigenstates via the control. This 
happens, for instance, in a planar rotating molecule controlled by only one control [BCCS12| 
Section 8] which is not (approximately) controllable. 

The use of more than one control opens new controllability horizons. 

Controllability results with more than one input have been obtained for specific sys- 
tems |EP09t IBBRIO] and some general approximate controllability results between eigen- 
functions have been proved via adiabatic methods |AB05| IBCMS12] . The first multi- input 
result via Lie-algebraic methods is given in |BCCS12| Section 8] where the problem of the 
spectral degeneracies in the planar rotating molecule has been overcome associating with ev- 
ery 1-dimensional slice of the set of admissible controls an invariant subspace of the state 
space l-L on which the single-input controllability result applies. Anyhow, such a technique 
does not apply for more general cases. In the case of a rotating rigid symmetric 3D molecule 
this application of this method is obstructed by the fact that eigenspaces may have arbitrarily 
large dimension, since at every higher energy level new degeneracies appear. In some sense 
the strategy of |BCCS12| Section 8] does not fully exploit the potentialities of the geometric 
method based on the controllability of the Galerkin approximations. 

In this paper, we present a sufficient condition for controllability of the discrete-spectrum 
bilinear Schrodinger equation which can be applied to cases in which the spectrum of the 
internal Hamiltonian Hq is very degenerate. The results fully exploit the use of more than 
one control and extend to simultaneous controllability, approximate controllability in if*, and 
stalking. Proving that a system is a stalker (i.e. it permits to track in modulus any given 
trajectory; for precise definitions see Section [2]) is a crucial issue when describing systems 
containing dissipative levels (and the dissipation is not taken into account in the mathematical 
model). In this case, a strategy is to keep the population of the dissipative levels as low as 
possible during the transitions in order to minimize the effects of the dissipation (see for 
instance the STIRAP model |CH90l IVHBBOll IBCG+02| 1 

The result presented in this paper is more general than those present in the literature even 
in the single-input case. For instance it applies to the Laplace-Dirichlet operator on a compact 
interval of M with a control term of the type u{t)x. Let us mention that in |BCCS12] . ap- 
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proximate simultaneous controllability of this model has been proved breaking the degeneracy 
between gaps among eigenstates through perturbation techniques. Here we prove the ap- 
proximate simultaneous controllability and stalking without perturbation arguments. In this 
framework, proving a controllability result without perturbation arguments allows to trans- 
late directly the constructive proof of the main result in an algorithm that provides explicit 
expressions of controls |CBCS11) . 



1.1 Brief description of the general results 

The main result of the paper is a sufficient condition for approximate simultaneous controlla- 
bility which we call the Lie-Galerkin Control Condition (see Definition 12. 5p . 

Roughly speaking, both the sufficient condition proposed in |BCCS12] and the one pre- 



sented here are based on the idea of driving the system with control laws in resonance with 
spectral gaps of the internal Hamiltonian Hq. However, whereas in |BCCS12] the only actions 
on the system obtained by resonance which are exploited for the controllability are those cor- 
responding to elementary transitions between two eigenstates, no such a restriction is imposed 
in the Lie-Galerkin Control Condition (see Section [T^ . 

The Lie-Galerkin Control Condition ensures strong controllability properties for the Galerkin 
approximations. Indeed it provides controllability for a fixed Galerkin approximation while 
avoiding the transfer of population to higher energy levels. This allows estimates on the dif- 
ference between the dynamics of the finite-dimensional Galerkin approximation and the ones 
of the original infinite-dimensional system. The Lie-Galerkin Control Condition also ensures 
a bound on the norm of the control achieving controllability which is uniform with re- 
spect to the prescribed tolerance, when the required transfer is between finite combinations 
of eigenstates. 

Under the Lie-Galerkin Stalking Condition, a slight modification of the Lie-Galerkin Con- 
trol Condition, we can prove that any trajectory can be tracked in modulus (see Theorem l2.8p . 

The Lie-Galerkin Control Condition under the additional assumption that the system is 
s- weakly coupled (see Definition 12. lip as introduced in |BCC] . allows to conclude that the 
system is controllable in H^^'^ (see Theorem 12. 12p . 



1.2 Application to the quantum angular momentum 

Rotational molecular dynamics is one of the most important examples of quantum systems 
with an infinite-dimensional Hilbert space and a discrete spectrum. Molecular orientation and 
alignment are well-established topics in the quantum control of molecular dynamics both from 
the experimental and the theoretical point of view (see [SS031 lSKA+04[ [SH06] and references 
therein). For linear molecules driven by linearly polarized laser fields in gas phase, alignment 
means an increased probability direction along the polarization axis whereas orientation re- 
quires in addition the same (or opposite) direction as the polarization vector. A large amount 
of numerical simulations have been done in this domain but the mathematical part is not 
yet fully understood. From this perspective, the controllability problem is a necessary step 
towards comprehension. 

We focus in this paper on the control by external fields of the rotation of a rigid linear 
molecule in R'^. This control problem corresponds to the control of the Schrodinger equation 
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on the unit sphere 5*^. We show that the system driven by three fields along the three axes is 
approximately controllable for arbitrarily small controls. 

Up to normalization of physical quantities (in particular, in units such that h = 1), the 
dynamics are governed by the equation 

^ dip{0,ip, ) _ ^\ _|_ (uJf') sin 6* COS + U2{t) sin 6* siny^ + u^it) cos 0)t/j(6, (p, t) (2) 

at 

where 6, if are the spherical coordinates, which are related to the Euclidean coordinates 
through the identities 



while A is the Laplace-Beltrami operator on the sphere S"^ (called in this context the angular 
momentum operator), i.e., 



The wavefunction ip evolves in the unit sphere S oil-i = L'^{S'^, C). 

As a consequence of the general multi- input result presented in Section [2] we have that ([2]) 
is approximately controllable with arbitrarily small controls. A stronger statement, including 
simultaneous controllability in and stalking, is given in Section [31 

Theorem 1.1. For every ip^ , tp'^ belonging to S and every 6 > 0, there exist T > and 
u e L°°([0,T], [0,(5]^), such that the solution ?/;(■) of equation corresponding to the control 
u and with initial condition ipQ, satisfies Wip'^ — 'ip{T)\\ < e. 

There are two main difficulties preventing the application to this system results previously 
in the literature. Firstly, we deal here with several control parameters, while those general 
results were specifically conceived for the single-input case. Notice that, because of symmetry 
obstructions, equation ([2]) is not controllable with only two of the three controls Mi, M2, M3. 
Secondly, the general theory developed in |CMSB09l[NerT0llBCCS"T2] is based on nonresonance 
conditions on the spectrum of the internal Hamiltonian. The Laplace-Beltrami operator on 
S*^, however, has a severely degenerate spectrum, since the £-th eigenvalue —ii{i + 1) has 
multiplicity 2i + 1. In |CMSB09] we proposed a perturbation technique in order to overcome 
resonance relations in the spectrum of the drift. This technique was apphed in |BCM+n9| to 
the case of the orientation of a molecule confined in a plane driven by one control. The planar 
case is already technically challenging and a generalization of the same technique to the case 
of three controls in 3D seems hard to achieve. 

1.3 Structure of the paper 

The structure of the paper is the following: in the next section we present the general multi- 
input abstract framework and the main abstract results. In Section Owe apply them to system 
([2]). The proofs of the abstract results are contained in Sections |U |5l and El 



X = sin 6 cos ip, y = sin 6 sin (p, z = cos 6, 



A 



1 d 
sin 9 89 




sin^ 9 dip"^ ' 
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2 Framework and main results 



Let p eN, 6 > 0, and U = Ui x ■ ■ ■ x Up with either Uj = [0, 6] or Uj = [-6, 6]. 

Definition 2.1. Let "H be an infinite-dimensional Hilbert space with scalar product (-, ■) and 
A, Bi, . . . , Bp be (possibly unbounded) skew-adjoint operators on "H, with domains D{A), 
D{Bi), . . . , D{Bp). Let us introduce the controlled equation 

^it) = {A + uiit)Bi + --- + Upit)Bp)^{t), u{t)eU. (3) 

We say that A satisfies (Al) if the following assumption is true: 

(Al) A has discrete spectrum with infinitely many distinct eigenvalues (possibly degenerate). 

Denote by $ a Hilbert basis {(pkjk&n of "H made of eigenvectors of A associated with the 
family of eigenvalues {i\k)km and let L be the set of finite linear combinations of eigenstates, 
that is, 

£ = y span{(/)i, . . . ,(/)fc}. 

fcGN 

We say that (A, Si, ... , Bp^ U, $) satisfies (A) if A satisfies (Al) and the following assump- 
tions hold: 

(A2) 0fc e D{Bj) for every G N, j = 1, . . . ,p; 

(A3) A + UiBi + ■ ■ ■ + UpBp : £ — 7- H is essentially skew-adjoint for every u eU . 

If (A, . . . , Bp, U, $) satisfies (A) then, for every {ui, . . . , Up) E U, A + uiBi -|- ■ ■ ■ -|- 
UpBp generates a subgroup e*(^+"i^i-i \-upBj,) q£ ^-j^g group of unitary operators U('H). It is 
therefore possible to define the propagator at time T of system ([2]) associated with a p- 
uple of piecewise constant controls u{-) = (mi(-), . . . , Up{-)) by composition of fiows of the type 
^t{A+uiBiA hupBp)^ moreover, . . . , Bp are bounded operators then the definition can be 
extended by continuity to every L°° control law (see |BMS82| Theorem 2.5]). 

Definition 2.2. Let {A, . . . , Bp, U, $) satisfy (A). We say that ([3]) is approximately con- 
trollable if for every ipo,ipi in the unit sphere of "H and every e > there exists a piecewise 
constant control function u : [0,T] ^ U such that Wipi — T^{iIjq)\\ < e. 

Definition 2.3. Let {A, . . . , Bp, U, $) satisfy (A). We say that ([3]) is approximately simul- 
taneously controllable if for every r in N, ipi, . . . ,ipr "H, T in U('H), and e > there exists 
a piecewise constant control u : [0,T] — )■ f/ such that 



< e, k = 1, . . . ,r. 



If, moreover, for every ipi, . . . E C and T G U('H) such that Tipi, . . . , Tip^ G C, there exists 
K > (not depending on e) such that u can be chosen to satisfy, in addition, \\u\\li < K, we 
say that ([3]) is -bounded approximately simultaneously controllable. 
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This last definition of controllability with a priori bound on the L^-norm of the control 
achieving controllability has been observed in preceding works |BCCS12| IChal2| . It implies 
a stronger controllability property as shown in Section 12.41 

Due to presence of the internal Hamiltonian and the boundedness of the controls, it is not 
possible in general to track, with arbitrarily precision, an unfeasible curve in S. We introduce, 
then, the notion of stalker, that is a system for which it is possible to track any given curve up 
to phases (both for a single initial condition and in the spirit of simultaneous control). This 
definition makes sense from the physical point of view, since tracking up to phases means 
imposing the population of all energy levels of Hq along the evolution. 

The identification up to phases of elements of "H in the basis $ = (0a:)a:gn can be introduced 
through the projection 



fcGN 



Definition 2.4. Let {A, Bi, . . . , Bp, U, $) satisfy (A). We say that ([3]) is a stalker if for every 
r in N, ipi,...,ipr in "H, T : [0,T] — )■ U('H) continuous, with Tq = Id^, and e > there 
exist an invertible continuous function r : [0,T] — >■ [0,Tt-] and a piecewise constant control 
u : [0, Tt] — )■ t/ such that 



< e. 



1, 



for every t G [0, T^]. 



2.1 Notation 

For every n in N, define the orthogonal projection 

n 
k=l 

Given a linear operator Q on H we identify the linear operator i^nQT^n preserving span{0i, . . . , 
with its n X 77, complex matrix representation with respect to the basis {(pi, . . . , We define 

= vr^Avr, and 5^"^ = vr^^.vr^, 

for every j = 1, . . . ,p. 

Let us introduce the set S„ of spectral gaps associated with the n-dimensional Galerkin 
approximation as 

= {|A/ - Afcl \ l,k = 1, . . . ,n}. 
For every a > 0, every m G N, and every m x m matrix M, let 

The nxn matrix Scr{Bj^^), j = 1, . . . ,p, corresponds to the "activation" of the spectral gap a: 
it reflects the action of the convexification procedure detailed in the following sections, which 
annihilates all the matrix elements {Bj^^)i f^ such that |A; — A^l a. 
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Define 



>, j) e S„ X {1, . . . ,p} I 3M e u(ri) s.t. £„{B 



M 








* 



for every N > nj . 

(4) 

Tlie matrices Sa{Bj"'^) for (a, j) G H„ correspond to "compatible dynamics" for the ?7,-dimensional 
Galerkin approximation {compatible, tliat is, with higher dimensional Galerkin approxima- 
tions) . 

2.2 Controllability results 

Let 

V° = U {^.(Sf ^) I (cr, j) e H„ and j is such that (0, j) G S„} 

U {^.(5]'^^) I G S„,a ^ 0, f/, = [-6,6]} . 

Definition 2.5. Let (A, Bi, . . . , Bp, U, $) satisfy (A). We say that the Lie-Galerkin Control 
Condition holds if for every no there exists n > uq such that 

LieV°Dsu(n). (5) 

Theorem 2.6 (Abstract multi-input controllability result). Assume that (A) holds true. If 
the Lie-Calerkin Control Condition holds then the system 

X = {A + UiBi + ■ ■ ■ + UpBp)x, u ^ U, 

is -bounded approximately simultaneously controllable. 

2.3 Stalking results 

For every ^ G S*^ C C, consider the matrix operator such that 

^Mj- fc if Xj < Afc 

(Mm,,, = {0 if A, = Xk 

^CMj^k if Xj > Afc. 



Let 



Vn = {7^(^.(5]"^)) I {a,j) G S„,a ^ 0,e G 5^} . 

Notice that V„, C 5u{n). 

Definition 2.7. Let {A, Bi, . . . , Bp, U, $) satisfy (A). We say that the Lie-Calerkin Stalking 
Condition holds if for every rio G N there exists n > uq such that 

LieVn = su{n). (6) 
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Theorem 2.8 (Abstract multi-input tracking result). Let Uj = [S,S] for some S > and 
every j = 1, . . . ,p. Assume that (A) holds true. If the Lie-Galerkin Stalking Condition holds 
then the system 

X = {A + uiBi + ■ ■ ■ + UpBp)x, u E U, 

is a stalker. 

Remark 2.9. If Uj = [—6, 6] for every j = 1, . . . ,p, then the Lie-Galerkin Stalking Condition 
implies the Lie-Galerkin Control Condition, as it follows from the relation 

2.4 Controllability in higher norms 

We define for s > 0, 

nGN 

for every belonging to 

D{\An = Len I J2 iAnri(^,0n)r < +ool . 

L nSN J 

For every G -D(|y4|**) we can define the |y4|*-norm (or simply s-norm) of tp by \\'ip\\s = 

If A is the Laplace-Dirichlet operator on some bounded domain of M" then the s-norm is 

equivalent to the if^*-norm on Z}(|yl|'*). 

Definition 2.10. Let {A, Bi, . . . , Bp, U, $) satisfy Assumption (A) and let s > 0. System ([3]) 
is approximately simultaneously controllable (respectively approximately controllable) for the 
s-norm if for every e > 0, r G N (respectively r = 1), ijji, . . . ,ijjr in -D(|A|^), and T G U('H) 
such that Ttpi, . . . , Ttpr G -D(|^|*) there exists a piecewise constant function : [OjT^] — )■ M 
such that 

for every j = 1, . . . , r. 

We say that [A, Bi, . . . , Bp, U, $) satisfies (A') if it satisfies (A) and the following additional 
assumptions hold: 

(A4) the operator i{A + UiBi + ■ ■ ■ ^- UpBp) is bounded from below for every u eW; 
(A5) the sequence {\k)km is non-increasing and unbounded. 

Definition 2.11. Let {A,Bi, . . . ,Bp,U,(^) satisfy Assumption (A') and let s > 0. Then 
{A, Bi,..., Bp) is s -weakly- coupled if D{\A + uiBi + ■■■ + UpBpY''^) = DdAl""/^) for every 
u eMP and there exists C such that 

m\A\^^,Bi^)\<C\{\A\^ij,^)\ 

for every / = 1, . . . ,p, ip G iI'(|A|'*). 
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The following result is a consequence of |BCCt Proposition 2] and can be obtained by 
adapting the arguments proposed in |BCCt Proposition 5]. We provide its proof in Section O 

Theorem 2.12. Let {A, Bi, . . . , Bp, U, $) satisfy Assumption (A') and (A, Bi, . . . , Bp) be s- 

weakly coupled for some s > 0. If is -bounded approximately simultaneously controllable 
then it is approximately simultaneously controllable for the s/2-norm. 

As a direct consequence we have the following result generalizing |BCCt Proposition 5]. 

Corollary 2.13. Let {A, Bi, . . . , Bp, U, $) satisfy Assumption (A') and {A, Bi, . . . , Bp) be s- 
weakly coupled for some s > 0. // the Lie-Galerkin Control Condition holds then system ^ 
is approximately simultaneously controllable for the s/2-norm. 

2.5 Example: the infinite potential well 

We present the case of a particle confined in the interval (—1/2,1/2) as a toy model to 
compare the result in |BCCS12] and Theorem 12.81 on a single-input system. The model has 
been extensively studied by several authors in the last decade and it has been the first quantum 
system for which a positive controllability result has been obtained (see |BC06) ). In |BCCS12 
an approximate simultaneous controllability has been obtained with geometric methods and 
using perturbations techniques. Indeed this model presents several resonances preventing the 
direct application of the results in |BCCS12] . 
The Schrodinger equation writes 

dip 1 d^ip 

with the boundary conditions ?/'(— 1/2, t) = tp{l/2, t) = for every t G M. We consider controls 
u{-) piecewise constant with values in U = [—S,6] for some S > 0. 

In this case "H = ((—1/2, 1/2), C) endowed with the Hermitian product {ipi,ip2) = 
J^y^ijji{x)il>2{x)dx. The operators A and B = Bi are defined by Aip = for every ip in 

D{A) = {H2 n Hi) ((-1/2, 1/2), C) and Bip = ixip. A complete set of eigenf unctions of A is 

, , , I v2cos(A;7rx) when k is odd, 
I v2sm(fc7rxj when k is even, 

,2 2 

associated with the eigenvalues i\k = —i—^, k eN. Notice that 

{<P„B<Pk)^0 

if and only if j + /c is odd. In particular V„ C su{n) for every n. 

We prove by induction on n that LieV^ = su{n), and hence that the Lie-Galerkin Stalking 
Condition is fulfilled. Notice that the matrices E2k-i{B^^^), for k < N, have only zero elements 
in the positions {j, I) for j < N and / > + 1, since 

/2 _ j2 > (iv + 1)2 _ = 2A^ + 1 > 2A; - 1. 

Hence S2k^i{B^^^) e Vn ior k = 1, . . . , N. We prove the claim by showing that 

Lie ({£2fc.i(5(")) \ k = l,...,n}) = su{n). (8) 
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For 77, = 2, 63(3^'^^) = y g^J generates su(2) because 612 7^ 0. 

Now assume that 

Lie ({£2fc-i(5^"-'^) \ k = l,...,n-l}) = su{n - 1), 

and let us prove ([8]). The matrices £2k-i{B^'^^) are in su(n) for every k = l,...,n — 1 
and generate the subalgebra of matrices in su(n) with zero elements in the nth row and 
nth column. In particular there exists M G Lie ({£^2fc-i : A; = 1, . . . , n - 1}) such that 

M + £^2n-i(-B has only two nonzero elements in the positions (n, — 1, n) and (n, n — 1). So 

Lie({<£:2fc-i(5("^)|A; = l,...,n}) 

D Lie ({£2fe-i(5(")) I A; = 1, . . . , n - 1} U {M + ^2n-i(5^"^)}) = su(n). 

Therefore, thanks to Theorems 12.61 and 12. 8[ system ([7]) is approximately simultaneously con- 
trollable and a stalker. 



3 The 3D molecule 

Let us go back to the system presented in the introduction for the orientation of a linear 
molecule, that is, 

ihip = —Aip + [ui cos 6 + U2 cos ip sin 6 + U3 sin ip sin 6)1^, (9) 

where ij{t) eU = L\S^,C)- 

A basis of eigenvectors of the Laplace-Beltrami operator A is given by the spherical har- 
monics Y^{6,ip), which satisfy 

AYr{e,ip) = -i{i + i)Yr{e,ip). 

The spectrum of A = iA is + 1) | ^ G N}. Each eigenvalue + 1) is of finite 

multiplicity 2i + 1. Therefore A satisfies Assumptions (Al) and (A5). Using the notations of 
the preceding sections we set Bi, B2, to be the multiplication operators by 
— i sine/) sin ^, —icos6 respectively. Being Bi, B2, B3 bounded, conditions (A2), (A3), and 
(A4) hold. Hence (A') is satisfied. Moreover, as proved in |BCCl Proposition 8], (Q is s- 
weakly coupled for every s > 0. The main goal of this section is to prove that system 
satisfies the Lie-Galerkin Stalking Condition. As a consequence we obtain the following result, 
whose corollary is Theorem ILII 

Theorem 3.1. System is: 

(i) -bounded approximately simultaneously controllable, 

(ii) approximately simultaneously controllable in for every s > 0, 

[iii) a stalker. 
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Using classical identities for Legendre polynomials and trigonometric relations one can 
prove that 

for every j = 1, 2, 3, and m, m' G {—I — 1, ...,£ + 1}. 
Moreover 

(yr,W') = o 

with \t - f I > 2 for every m G {-£ - 1, . . . , £ + 1}, m' G {-f - 1, . . . , f + 1}, j = 1, 2, 3. In 
order to prove that the Lie-Galerkin Stalking Condition is satisfied, we choose a reordering 
{(j)k)k<m of the spherical harmonics in such a way that 

I A: = l,...,4£ + 4} = {y,-^...,F/,r-^^^•••,>^m}, 
and we are left to prove that 

LieV4m = su(4£ + 4). 

The characteristics spectral gap of the space l-Li is (£ + 1)(^ + 2) — l{l + 1) = 2{l + 1). In 
particular (2(£ + 1), 1), (2(£ + 1), 2), and (2(£ + 1), 3) are in ^4^+4. 

3.1 Matrix representations 

Denote by the set of integer pairs {(r, m) | r = ^, £ + 1, m = — j, . . . , j}. Consider the 
lexicographic ordering q : {1, . . . , 4£ + 4} — >■ JT^. For j, = 1, . . . , 4£ + 4, let e^^fc be the (4£ + 4)- 
square matrix whose entries are all zero, but the one at line j and column k which is equal to 
1. Define 

Bj,k ^i,k ^k,ji Pj,k ^^j,k ~l~ i^kj) ^j,k ^(^k,k- 

By a slight abuse of language, also set eQ{j),e{k) = Gj,k- The analogous identification can be 
used to define -Dg(j),p(fc). Note that 

<^i(-E'(^,m),(f+l,n)) = —F'{e,m),{e+l,n), and Ji{F(^l,m),(£+l,n)) = -E'(f,m),(£+l,n)- (10) 

Thanks to this notation we can conveniently represent the matrices corresponding to the 
controlled vector field (projected on He)- A computation shows that the control potentials in 
the X and y directions, — i cos </? sin ^ and —isimpsinO respectively, projected on Tie, have the 
matrix representations 

e 

m=-e 

^ ^ = {(le,mE(l^m),{e+l,m-l) + (le~m.E{e,m),{e+i,m+i)), 

m=—£ 

where 

/(£-m + 2)(£-m + 1) 
^'"^ ~ V 4(2£ + l)(2£ + 3) ■ 
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Similarly, we associate with the control potential in the z direction, —% cos Q the matrix 
representation 

I 

r(4£+4) _ „ p 

-D3 — 2^ P^,m-r(£,m),(^+l,m), 

m=—l 

with 

(^ + 1)2 -m2 
_ (2^+1) (2^ + 3)' 

3.2 Useful bracket relations 

From the identity 

we get the relations 
and 

[-E'i.fc, -^i,fe] = SZ^j^fc, and [Fj^k, Dj^k] = '^Ej^k- (12) 

The relations above can be interpreted following a "triangle rule": the bracket between an 
operator coupling the states j and k and an operator coupling the states k and n couples the 
states j and n. On the other hand, the bracket is zero if two operators couple no common 
states, that is, 

fc, Z,, = if {j. A;} n {f, k'} = 0, (13) 

with Y,Z e {E,F,D}. 



3.3 Controllability in 5u{4£ + 4) 



Lemma 3.2. The Lie algebra L generated by i?f i?f J.(5f ^^(^f 

Ji{Bf^^^^) is equal to su(4£ + 4). 

Proof. The first step of the proof consists in showing that the Lie algebra L contains the 
elementary matrices 

E[i,k),{i+i,k+3) for A; = -£,... , j = -1, 0, 1. (14) 

With a slight abuse of notation and for the sake of readability, let us write Bj = B^^^^^\ 
j = 1, 2, 3. Let us also write ado, /3 for [a, /3] and ad;^"*"^ j3 for [a, ad;^ f5]. 



Notice that 



■JiiBs) — Pl,mJi{F(£^m),(e+l,m)) — Pe,mE(£^m),(e+l,m) ■ 

£=—m £=—m 
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By induction on j > and using the bracket relations fll2p , we have 

ad| MBs) = [Bs, [53, ad^-2 MB^)]] 

m=-i 

By invertibihty of the Vandermonde matrix and since pi^m 7^ Pe,n for every n ^ m, — m, it 
foUows 

E(e-m),(e+i-m) + E(^i^rn),{e+i,m) G L for m = 0, . . . , (15) 
In particular -E{£,o),{£+i,o) ^ L. The double bracket of 

B2- Ji{Bi) _ ^ T n«^ 

m=—£ 

with -E'(£,o),(€+i,o) is easily computed using (ITT]) and (fT3|) and gives 

^ ^ " '^^(A0),(A1)5 -C/(f,0),{£+l,0) 



1^1^ Q'f,m-E'(f,m),(£+l,m-l)5 -E'(£,0),(£+1,0) 5 -E'(€,0),(£+1,0) = [-E'(£,0),(A1) ! -E', 



m=-£ 



~ (ltfi[E{i+l~l),{t+lfi)-, -£'{f,0),(f+l,0)] 
= 9£,0-E'(€,0),(f+l,-l) + 5'£,i-^(^,i),(^+i,o) ^ 

Define Qq = g£,o-E^(£,o),(£+i -i) + g£,i-^(^,i),(^+i,o) and, similarly, = qi^.jnE(e~m.),{e+i~m~i) + 
qe,m+iE(i^rn+i),[£+i,'m) foY < iTi < i, Qe = g£-£-E(£ _£),(£+! . In particular B2 - Ji{Bi) = 

Using again (ITT]) and f lT3|) . we have 



m=A: 



for /c = 1, . . . , By recurrence on k and because of (I15p . it follows that Qk & L for /c = 0, . . . , 
Now, since Qi/qi-i = _£)_(£+i is in L, then 

which, in turns, implies that 

Iterating the argument, E(^i^rn),i£+i,m) and -E(£^m)^(<>+i^„_i) are in L for every m = -£,... 
Developing the same argument as above replacing f fT6|) by 



m=—£ 



we have that also -E'(£,m),(^+i,m+i) is in L for every m = —£, proving (fT^ . It then follows 

from ( II ip that each Ej^^ is in L. 
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If now we replace fll6p by 



Bi + JiiB2) _ p j- 

^ — ~ Qe,m-r {e,m),{e+l,rn-l) ^ J^, 



m=-i 

or 



Bi - Ji{B2) _ p J 



m=-l 



we obtain from the arguments above that F(i^rn),{i+i,m~i) and F(i^rn),{i+i,m+i) are in L for every 
m = -i, . . . ,i. The relations ([II]) and ([l2]) allow then to conclude that L = su(4£ + 4). □ 



4 Proof of Theorem 2.6 



4.1 Time-reparametrization 

Up to replacing each 5^ by 6Bj, we can assume that 5 = 1. 

For every piecewise constant function z such that z{t) > 1 for every t, we consider the 
time-reparametrization 

^(t) = izit)A + + ■ ■ ■ + Upit)z{t)Bp)ij{t) (17) 

of system Each Uj{t)z{t) belongs to the time- varying set z{t)Uj. 

If Ml, . . . , Mp are control laws in (fT7|) then the corresponding controls in ([3]) are their time- 
reparametrizations Uj{s) = Uj(t{s)) with = z{T)dT, j = 1, . . . ,p. By restricting the 
range of available controls and setting Vj(t) = Uj{t)z{t), we can focus our attention to trajec- 
tories of 

^(t) = izit)A+v,{t)B^ + ---+Vpit)Bp)^Pit), z{t) > 1, v{t) = {v^{t),...,vp{t)) G U. (18) 

Each solution of (fT8|) with z and f piecewise constant is the time-reparametrization of a 
solution of ([3]) with piecewise constant controls (but the converse is not necessarily true, 
since we restricted the set of admissible controls). Hence, the approximate simultaneous 
controllability of ( I18p implies the approximate simultaneous controllability of Moreover 



\u 



\uj{T)\dT = / \uj{t)\z{t)dt = / \vj{t)\dt<T, 



for j = l,...,p. The last inequality holds since either Uj = [0,1] or Uj = [—1,1]. Hence 
the approximate simultaneous controllability in C of (fTSj) with a bound on the controllability 
time uniform with respect to the tolerance implies, in fact, the L^-bounded approximate 
simultaneous controllability of 
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4.2 Interaction framework 

Given a solution of f|T8|) with controls z{-), fi(-), . . . , Vp{-) and a piecewise constant func- 
tion a(-) with values in {0, 1}, let us define 

^i't) = / i^i^) ~ 0!{s))ds 
Jo 

and 

In particular 

|(0fc,y(t))| = |(0,,^(t))|, A;eN, (19) 
for every t. For f i, . . . , G M set Q{uj, Vi, . . . , Vp) = e~'^'^{viBi + ■ ■ ■ + VpBp)e^^. Note that 

e(a;, vi,..., Vp)jk = {(pk, v^, . . . , Vp)(f)j) = e*{^'=-^^)- {v^{Bi)^k + ■■■ + Vp{Bp)^k) , (20) 

and that y{-) satisfies 

yit) = iait)A + Q{uit),viit), . . . , Vpit)))y{t), a e {0,1}, v e U,u + a > 1. (21) 

Conversely, each solution of (12T]) with a G {0, 1} and v & U piecewise constant and u contin- 
uous and piecewise affine, with u + a = z > 1 almost everywhere, is, up to a time-dependent 
change of coordinates preserving the modulus of each component with respect to the basis $, 
a solution of (llSp with u piecewise constant. In particular, each solution of 

yit) = {a{t)A + eiLo{t),v,{t), . . . , Vp{t)))y{t), a G {0, l},veU,io> 1, (22) 

with a, V piecewise constant and u continuous and piecewise affine is, up to a time-dependent 
change of coordinates preserving the modulus of each component, a solution of (ITS]) with u 
piecewise constant (but the converse is not necessarily true). 

Proposition 4.1. Approximate simultaneous controllability of fl22]) implies approximate si- 
multaneous controllability of ([3]). //, moreover, approximate simultaneous controllability in 
C = [J^gj^ span{0i, . . . , of (122|) is achieved with a uniform bound on time then (j3]) is 
-bounded approximate simultaneous controllable. 

Proof. The strategy of the proof follows the idea of the proof of |BCCS12| Proposition 6.1]. 

It follows from (fT9|) that approximate simultaneous controllability of (122|) implies approx- 
imate simultaneous controllability of (ITS]) in modulus. 

Moreover, because of the unitarity of the evolution, the approximate simultaneous con- 
trollability of (122|) is equivalent to the approximate simultaneous controllability of the system 

yit) = -{a{t)A + Q{uit),vi{t),...,Vp{t)))y{t), aG{0,l}, v e U, 

which implies approximate simultaneous controllability in modulus of the time- reversed of (fTSjl . 

Take r orthonormal initial conditions iIjq,...,iPq and r orthonormal final conditions ipl, . . . 
Since the the spectrum of A is infinite by Assumption (Al) we can apply |BCCS12| Lemma 6.3] 
so that for every tolerance 77 > there exist ki, . . . ,kr such that 

C = {e*^0fc, + ■ ■ ■ + e*^0fc, I t G M}. 
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is ?7-dense in the torus 

r = {e^i^0fc, + ■ ■ ■ + e^'-^0fc, I ^1, . . . , G M}. 

By approximate simultaneous controllability in modulus of (fT8|) it follows that there exists 
an admissible control (2;, v) steering simultaneously each ip^, for j = 1, . . . , r, ?7-close to e^^^(j)k^ 
for some 61, . . . ,6r & ^■ 

Similarly, by approximate simultaneous controllability in modulus of the time-reversed 
of f lTSj) there exists an admissible control [z, v) steering system (fT8|) simultaneously, for some 
9i, . . . , 9r & from e^^'^^fcj, . . . , e^''^(j)k^ to an ?7-neighborhood of , . . . , 

Finally the concatenation of the control {z,v), a control constantly equal to (1,0) on a 
time interval of suitable length, and [z, v) steers system ffTSjl simultaneously from ipl, . . . 
to a 3?7-neighborhood oi ipl, . . . 

According to the conclusion of Section 14. 1^ the approximate simultaneous controllability 
of (1T8|) implies approximate simultaneous controllability of (|3]). □ 

4.3 Galerkin approximation 

Definition 4.2. Let N eN. The Galerkin approximation of (122|) of order is the system 

X = {aA^^^ + e^^\uj,vi,...,vp))x, xen, (23) 

where Q^^^co, Vi, . . . , Vp) = 7riyQ{uj, Vi, . . . , Vp)txn. The controls v are piecewise constant with 
values in f/, while u is continuous and piecewise affine, with Cj >1 almost everywhere. 

In the following section we recall a convexification result whose role is to identify the 
matrices that can be obtained by convexification of matrices of the form G*^^-'(a;, f 1, . . . ,fp). 
Recall that the elements of 6^-^)(ci;, Ui, . . . , fp) are described by (120]) . 

4.4 Convexification 

The following technical result has been proved in |BCCS12 



Lemma 4.3. Let k he a positive integer and 71, . . . , 7^ G M \ {0} he such that |7i| ^ |7j| for 
j = 2, . . . , n. Let 

ip{t) = (e^*^S...,e^*T-). 

Then, for every tq G M, we have 



convv9([ro, 00)) 3 CDS^ x {(0, . . . , 0)} , 

where 

00 

^ = n^°^©>0- (24) 

k=2 

Moreover, for every R > and ^ E there exists a sequence (rfc)^-^ such that Ti > Tq, 
Tk+i -Tk> R, and 

1 ^ 

hm ^^(r,) = (0)^,0,. ..,0). 

k=l 
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4.5 Choice of the order of the Galerkin approximation 

In order to prove approximate simultaneous controllability, we should take r in N, ■j/'i, . . . , V'r in 
"H, T in U('H), and e > and prove the existence of a piecewise constant control u : [0, T] — )■ f/ 
such that 

ftpk - Tlj^tpk <e, k = l,...,r. 
Notice that for uq large enough there exists g G SU{nQ) such that 

I ^"^fe) - (7r„o0i, grrnoipk) \ < e (25) 

for every 1 <k < r and j G N. This simple fact suggests to prove approximate simultaneous 
controllability by studying the controllability of the lift of (122|) in the Lie group SU{no). 

4.6 Control in SU{n) 

Let n > no be chosen, in accord with the statement of Theorem \2.6\ such that hypothesis ([5]) 
holds true. Define the set of matrices 

>V. = {A(")}u{fo(i?f^)l(o,j)e2„} 

U 1^0(5]"^) + CD^^(5]"^) I ((T,j) e En and a,j are such that (0, j) G S„,(T ^ o| 

U |cd£,(5]"^) I (a, j) G S„, or ^ 0, and f/,- = [-1, 1]} , 

where H„ and W are defined as in and flM|) . respectively. (Recall that by rescaling we are 
assuming 6=1.) 

Notice that Lie(W„) = Lie(V°). 

Consider the auxiliary control system 

X = M{t)x, M{t) G >V„, (26) 

where M plays the role of control. It follows from ([5]) and standard controllability results on 
compact Lie groups (see |JS72) ) that for every g G SU{n) there exists a piecewise constant 
function M : [0, T] W„ such that 



exp / M{s) ds = g, 
Jo 



where the chronological notation exp Vs ds is used for the flow from time to t of the 
time- varying equation q = Vs{q) , q E C"' (see |AS04] ). 

4.7 System reduction by convexification 

Let n be fixed as in the previous section. For every > n let 

W.,^ = {A(^)}u{^o(fir)l(o,j)eH„} 

U ^So{Bf^) + COS^{Bf^) I {a,j) G En and a,j are such that (0, j) G H„,a ^ o| 
U {wS^iB^) I (a, j) G S„, a ^ 0, and U, = [-1, 1]} . 
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Lemma 4.4. For every N > n and for every piecewise constant M : [0, T] — )■ yVn,N there exist 
a : [0,T] — )■ {0, 1}, v : [0,T] — t- U piecewise constant and a sequence {oJh{'))heN of continuous 
and piecewise affine functions from [0,T] to [0, oo) with ujh > 1 almost everywhere, such that 



uniformly with respect to t E [0, T] as h tends to infinity. 



M{s)ds 







Proof. Let N > n. Let us fix a and Vi, . . . ,Vp at each t E [0,T] as follows: if M(t) = A^^^ 
then a{t) = 1 and vi{t) = ■■■ = Vp{t) = 0; otherwise, if M{t) = So{bP), M{t) = £o{Bf^) + 



COScr{Bj^^), or M{t) = wScr{Bj" ') for some j, then take such a j minimal and set Vj{t) = 1 
and a{t) = Vk{t) = for A; 7^ j. 



AN) 



We are going to apply Lemma 14.31 for every interval on which M(-) is constant. Fix 
u}h{0) = for every h. Take an interval (to)^i) on which M(-) is constant and assume that 
U!h{to) has been computed. We next extend cUh on (to,ti). 

If a = 1 on (to, ^i) then take Uh^r) = cohito) + t — to for every r G (to, ti). 

Otherwise, let Vj = 1 on (to^^i) and assume for now that M{t) = £o{Bj^^) + w£rj{Bj'^^). 
Apply Lemma lOl with 71 = cr, {72, . . . , 7k} = ^n\ {o"}, = ^, R = T, and tq = Uh{to)- Then 
there exists a sequence {Tk)'^i such that ri > uJh{to), Tfc+i — Tk > T, and 



lim ^y^{e' 



3-^=7.) = (CO, o,...,0). 



fe=i 



In particular there exists K = K[h) such that 



fc=l 



< 



for every 1 < l,m < N. 

Consider the piecewise constant function Y : (to,^i) — defined as follows: set Sc 
ti + (^1 — to)a/K, a = 0, . . . , K, and let 



K 



Y{t) = ujh{to) + ^r„X[.c-i,sc)W • 



Following the smoothing procedure of |BCCS12| Proposition 5.5] one can construct a con- 
tinuous piecewise affine approximation Uh '■ [^Oj^i] — > M of F with > 1 almost everywhere 
such that 



(eW(u;,(.),i;i(^ 



,Vp{s)))ds- / M{s)ds 



ta 



ta 







(27) 



uniformly with respect to t G [to,ti] as h tends to infinity. 

The same argument can be carried out in the case in which M{t) = SoiBf^) by applying 
Lemma SSlwith 71 in (0, 00) \ T.^, {72, ■■■,1k} = ^n, ^ = R = T, and Tq = W/,(to)- 
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The final case to be considered is when M{t) = (x}Scr{Bj^^) with a 7^ 0, (c, j) G S„, and 
Uj = [-1, 1]. Notice that 

The argument above can be easily adapted to matrices M{t) of the type Vj{SQ{Bj^^) + 
G)J^{£^{bP))), with Vj e Uj, e e 51 (just not imposing ^ = 1 while applying Lemma [4.3p . 
and in particular to -{So{Bf^) + wJ.i{S^{Bf^))). 

It suffices then to introduce a sequence {M^')h^fq of piecewise constant functions with values 

in 

{v,{So{Bf^) + WJ^iS^B^^))) I V, G f/„ e G 5^} 

such that M^{s)ds converges uniformly for t G [to^^i] to J^^ M{s)ds as /i tends to infinity 
and to apply a diagonal procedure based on the approximation introduced above. □ 

As a consequence of the lemma above and thanks to |AS04t Lemma 8.2], we have 



exp 



(29) 



uniformly with respect to t G [0, T] as h tends to infinity. 



4.8 Control of the infinite-dimensional system 

Next proposition states that we can pass to the limit as tends to infinity without losing the 
controllability property fl2^ . Its proof is based on the special sparsity structure of the matrices 
in yVn,N, guaranteeing that the difference between the dynamics of the infinite-dimensional 
system and the dynamics of the Galerkin approximations is small. 

We introduce the following notation: given n G N and a bounded linear transformation L 
of "H, let Crop„(L) be the n xn matrix ((0^, L0fc))"^^^. We use the same symbol Crop„ to 
denote the similar cropping operation acting on the space of N x N matrices, with N > n. 

Proposition 4.5. Let n G N and M : [0,T] — Wn be piecewise constant. Then, for every 
5 > 0, there exist piecewise constant controls z : [0, T] — > [l,oo) and v : [0,T] U and a 
continuous piecewise affine function u with ci; > 1 such that the propagator \I' of (I22p satisfies 



exp f M{s) ds - Crop„(^ft) 
Jo 



for every t G [0, T] . 



Proof. Consider /i > to be fixed later. For every j G N the hypothesis that (pj belongs 
to D{Bi) implies that the sequence {{Bi)jk)k(zfq is in i"^ for every / = 1, . . . ,p. It is therefore 
possible to choose N > n such that \\{{Bi)jk)k>N\\e'2 < fJ' for every j = 1, . . . ,n and I = 1, . . . ,p. 

Let M be a piecewise constant function from [0,T] to yVn,N such that Crop„ M(t) = M{t) 
for every t in [0,T]. Because of the definition of S„ and of the classes Wn,N and W„ we have 
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exp / M{T)dT 



exp / M{T)dT 








* 



By Lemma r4.4[ for every > 0, there exist piecewise constant controls a : [0,T] — {0, 1}, 
V : [0, T] — )■ f/ and a continuous piecewise affine function u with uj >1 such that 



e^ / (a(r)A(^) + e(^)(w(r),t;i(r),...,?;p(r))) rfr- e"^ / M(r) rfr 



for every s,t in [0, T]. 

Consider the solution \E' of fl22|) associated with a, w and v. Set, for /c G N, 



gf)=Crop,vl/, 



Now 



(30) 



and the choice of is such that 

for every j = 1, . . . , n and k = 1, . . . ,N. Notice, moreover, that the norm of -RI^"* can be 
uniformly bounded by a positive constant C independent of rj (possibly depending on and 
hence on fi). 

By the variation formula and since Q[^^ = Crop„((5l^'') we have 

-t 



qJ") = Crop, 



e"^ / (a(r)A(^) + 0(^)(a;(r),i;i(r),...,^,(r))) cir + 



\ 



(^p r(«(r)AW + e(^)(a;(r),t;i(r),...,t;p(r))) dr) RfU^ 



so that 



CiopJ^-^p / {aiT)A^''^ + e^''\uiT),v,iT),...,VpiT))) dr 



< triC 

The norm of the matrix product 



(^e^p M{t) dr^ Crop„ {R^^) ds 



e^p f M(r) Crop„ (i?^^)) 

J s 



is equal to 



l|Crop„/?W||. 

The max norm of Crop„ i?!^^ is smaller than jj as it follows from fl30|) . Hence 

Ctop^ (^%- e^£(a(r)A(^) + e(^)(u;(r),t;i(r),...,t;,(r))) dr 
< T{r]C + V^fi). 
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The constant T{riC + ^/nf^) can be made arbitrarily small by choosing fi small with respect 
to n and T and then rj small with respect to C = C(/i) and T. □ 

Proof of Theorem \2.(A Let r in N, . . . , ■i/'r in "H, T in U('H), and e > 0. Let n be as in 
Section 14.61 and let g G SU{n) satify (125|) . Notice that if t/;!, . . . , t/^j.; 't('^i); • • • ? ^ ("^r) are in 
C then n can be taken independently of e. 

From Section 1121 there exists M : [0, T] — j- W„ such that 



exp 



M{s)ds = g. 



Proposition 14.51 ensures the existence of two piecewise constant functions z and v and of a 
continuous piecewise affine function u with uj > 1 almost everywhere such that the associated 
propagator \E' of (122|1 satisfies 



exp f M{s) ds - Crop„(^'T) 
Jo 



< e. 



If ipi, . . . , ijjr, T('?/'i), . . . , T('?/'r) are in £ then T is independent of e. By Proposition 14.11 sys- 
tem ^ is L^-bounded approximately simultaneously controllable. □ 



5 Proof of Theorem 12.8 

In order to prove Theorem 12.81 we adapt the proof of Theorem 12. 6[ The key point of the 
argument is the following: it has been proved in Proposition 14.51 that system (llSp can track 
every trajectory of (126|) . The idea is to replace f l26|) by a system which can track with arbitrary 
precision every trajectory in SU{n). The crucial property, beyond the Lie bracket generating 
condition, that the new version of fl26|) should satisfy in order to achieve this goal is that 
it is a driftless system (i.e., the time- reversal of each of its admissible trajectories is itself 
admissible) . 

The same time-reparameterization and time-dependent change of coordinates as in Sec- 
tion |4]T] allows to consider the tracking problem for system (122|) instead of system (|3]). As in 
the previous section we consider 5 to be renormalized to 1. 

We can then base our argument on the following analogue of Proposition 14.11 

Proposition 5.1. // (122|) is a stalker then ^ is a stalker as well. 

The following proposition allows to reduce a tracking problem in the space of unitary 
operators of % into a tracking problem in SU{n) for n large enough. Its proof can be found 
in |B(X:S12[ Proposition 5.7]. 

Proposition 5.2. Let T : [0, T] — )■ U('H) be a continuous curve. Take e > and m G N. 
Then for n > m sufficiently large there exists a continuous curve g : [0, T] — )■ SU{n) such 
that I {(pj, T(t)0fc) — {cj, g(t)ek)\ < e for every t in [0,T], 1 < k < m, and j = 1, . . . ,n, where 
Ci, . . . , e„ denotes the canonical basis o/M". 
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Let n be chosen as in Proposition 15.21 In accord with the Lie-Galerkin Stalking Condition, 
we can assume, without loss of generality, that Lie(Vn) = 5u(n). 

The roles played in Sections 14.61 and 14.71 by Wn and Wn,Af are now played by G)Vn and 
COVn.AT with 

In particular, we consider as auxiliary control system 

X = M{t)x, M{t) e CDVn, (31) 

M being the matrix-valued control parameter. It follows from the equality Lie(V.„) = su{n) 
and Rashevski-Chow's theorem that every trajectory on SU{n) can be tracked with arbitrarily 
precision (up to time-reparameterization) by a trajectory of (13T]) . 

The relation between the trajectories of (123|) and those of (13T]) (or, more precisely, x = 
M{t)x, M{t) G CDVn,7v), is described by the following lemma. 

Lemma 5.3. For every N > n and for every piecewise constant M : [0,T] — > wVn,N there 
exist a : [0,T] — ?■ {0,1}, v : [0,T] ^ U piecewise constant and a sequence {uh{-))hm of 
continuous and piecewise affine functions from [0, T] to [0, oo) with 0Jh> ^ almost everywhere, 
such that 



(a(s)A(^) + 0(^) {uh{s), v,{s),..., vp{s))) ds - f 

Jo 

uniformly with respect to t ^ [0, T] as h tends to infinity. 



M{s)ds 







Proof. The proof is almost identical to that of Lemma 14.41 in the case M{t) = {x)£rj{Bj 
The only difference is in replacing fl28|) by 

^j... ^iN),, (go(i^f ^) + g^MSAbP))) - {So{bP) + 0) J_g(g.(<^))) 

^Jd^a[Bj )) = . 

We then apply the same convexification argument. 

As in the previous section, the lemma above and |AS04t Lemma 8.2] imply that 



(7V)^ 



□ 



exp 



{a{s)A''^^ + e'-^\uJh{s),Vi{s),...,Vp{s))) ds- exp [ M{s)ds 

Jo 







uniformly with respect to t G [0,T] as h tends to infinity. 

In analogy with Section 14.81 we can conclude the proof of Theorem 12.81 thanks to the 
proposition below, which states that we can pass to the limit as tends to infinity without 
losing the tracking property of the finite-dimensional Galerkin approximations. Its proof is 
basically the same as that of Proposition 14.51 

Proposition 5.4. Let n G N and M : [0,T] — t- £DV„ be piecewise constant. Then, for every 
e > 0, there exist piecewise constant controls z : [0,T] — )■ [l,oo) and v : [0,T] — U and a 
continuous piecewise affine function u with a; > 1 almost everywhere such that the propagator 
^ of (122]) satisfies 



exp [ M{s) ds - Crop„(^t 
Jo 



for every t G [0, T]. 
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Remark 5.5. The hypothesis that each Uj contains in its interior can be relaxed. Indeed, 
up to reordering, let p' be such that Uj = [0,6] if j = 1, . . . ,p' and Uj = [—6, 6] for j > p'. 
Assume that, for every j G {1, . . . ,p'}, ii I ^ k are such that A/ = Afc, then {(pi, Bj(j)k) = 0. 
Assume, moreover, that the Lie-Galerkin Control Condition is satisfied with Vn replaced by 
the set of all matrices J^(£^o-(-B]"'*)) with (a, j) G H„, ^ G S^, a ^ 0, and either j > p' or the 
following holds: if /, k, I', /c' G {1, . . . , n} satisfy A; — A^ = A^/ — A^/ = cr and 

{(Pi,Bj(Pk) j^Oj^ {(P,,B,^k') 

then 

{4>i,Bj(f)i) - {(f)k,Bj(j)k) = {(f)i',Bj(f)i>) - {(j)k',Bj(f)k')- 

In this case the proof of Lemma [5 . 3 1 b ecomes more technically involved. The point is that, 
even if £^o(-Sj^'') cannot be eliminated by convexification, it is a diagonal matrix by hypothesis. 
Hence, it can be used to define a new interaction framework. The sequence {uh{-))h can then 
be constructed by following |BCCS12| Proposition 5.5]. 



6 Proof of Theorem 2.12 



First, let us prove L^-bounded approximate simultaneous controllability in s/2-norm for initial 
and final data in C. Namely we want to prove that, for r & N, tpi, . . . ,ipr & C, and T G U('H) 
with Ttpi, . . . , Ttpr G C, there exists K > such that the following holds: For every e > 
there exists a control u, with \\u\\li < K such that 

U3-^T^M/2<e j = l,...,r. (32) 

Let be such that ipi, . . . and Tipi, . . . , Tipr are in span{(/)i, . . . , (p^}- Note that on 
span{0i, . . . , 0Ar} we have 



N 



1/2 



s/2 



\hmk,^P)\' < (max{|Ai|, |A^|})^/2||V'||. (33) 



Since the system is L^-bounded approximately simultaneously controllable there exists A" > 
such that for every s > there exists a piecewise constant control u, with < K such 

that 

Wijj - TJi.ijj\\ < e j = l,...,r. 

Hence, by ([23D, we deduce L^-bounded approximate simultaneous controllability in s/2-norm 
in C. 

Now let ipi,...,ipr e D{\A\'/'^), and t G V{n) be such that t^i, . . . , t^^ G D{\A\'/'^). 
Let e > and consider ip^ , . . . , and ipl , . . . , ipl in C such that 

Wi'j - i'j\\s/2 < £ and - fi;j\\s/2 < e, 

for j = 1, . . . ,r. As proved above (see fl32|) ). there exist K, independent on e, and a piecewise 
constant control u with \\u\\li < K such that 



,1 'V^U 



r^^'^\\s/2<e J = l,...,r. 
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By |BCCt Proposition 2], since the system is s- weakly coupled, there exists a constant C 
depending only on s and A, Bi, . . . , Bp such that 

\\r^T^\\s/2 < C^||^||./2, 

for every e D{\A\'/^). Therefore 

< (C^ + 2)£, 

for j = 1, . . . , r. 

Remark 6.1. Using arguments similar to those of the proof of Theorem 12.61 and of Theo- 
rem 12.121 it is possible to prove a finer statement than Corollary 12.131 Indeed it is possible to 
prove that a system satisfying Assumptions (A'), the Lie-Galerkin Stalking Condition, and 
which is s-weakly coupled is a stalker for the s/2-norm. This is due to the fact that, actually, 
the Lie-Galerkin Stalking Conditionimplies stalking in C with a uniform bound on the 
norm of the control. 
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